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Abstract: We have searched for the decay Kg — > 77 in a sample of ~ 2 x 10^ 4> — > KsK^ 
decays collected at DA^>NE with an integrated luminosity of 1.9 fb~^. Ks are tagged 
by the Kl interaction in the calorimeter. Two prompt photons must also be detected. 
Kinematic constraints reduce the initial 6x10^ events to 2740 candidates, from which a 
signal of 711 it 35 events is extracted. By normalizing to the Ks — > decays counted 
in the same sample, we measure BR{Ks 77) = (2.26 ± 0.12stat i 0.06syst) x 10~^, in 
agreement with O(p^) Chiral Perturbation Theory predictions. 



Keywords: 



£+e Experiments. 



Contents 



|l]. Introduction 1 

m. The KLOE detector 2 

|3|. Search of Ks — ^ 77 with a pure Ks beam 2 

|3.1| Ks tagging and event preselection 2 



3.2 Kinematic fitting and event counting 4 

^ Branching ratio evaluation and systematics 7 

|5]. Conclusion 9 

|6|. Acknowledgements 10 



1. Introduction 

A precise measurement of tlie Ks^j'j partial width provides a test of Chiral Perturba- 
tion Theory (^PT). The Ks^2j decay amplitude has been evaluated at leading order in 
xPT, 0{p'^), providing an estimate to a few percent accuracy of branching ratio (BR): 
BR(i^'5'— >27)= 2.1 X 10~^ Q. Measurements of such BR have changed considerably with 
time m while improving in precision. The latest determination comes from NA48 Q, 
BR=(2.71 ±0.07) X 10"^ This resuh differs by about 30% from the £>(/) xPT estimate, 
possibly due to higher order corrections. 

We report in the following on a measurement based on a integrated luminosity J Cd t ~ 
1.9 fb"^ collected with the KLOE detector g at DA^>NE [§, the Prascati ^-factory 
DA$NE is an e^e~ collider operated at a center of mass energy, W, of ~ 1020 MeV, 
the mass of the 0-meson. Equal-energy positron and electron beams collide at an angle of 
(7r-0.025) radians, producing (/>-mesons nearly at rest. (/)-mesons decay 34% of the time into 
nearly collinear K^K^ pairs. Since J^'-''((j)) = 1 , the K^K^ pair is in an antisymmetric 
state and the two kaons are always a pure RgK^ pair. Detection of a i^T^-meson therefore 
guarantees the presence of a X^-meson of known momentum and direction. This proce- 
dure, called tagging, allows us to obtain a pure Kg beam. The data analyzed consists of 
some 2 billions KsKl pairs. 

* Corresponding author. Email address: matteo.martini@lnf.infn.it 
t Corresponding author. Email address: stefano.miscetti@lnf.infn.it 
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2. The KLOE detector 



The KLOE detector consists of a large cylindrical drift chamber, DC 0, of 4 m diameter 
and 3.3 m length operated with a low Z and density gas (helium-isoC4HiO), surrounded by 
a lead-scintillating fiber calorimeter, EMC [Q]. The chamber provides tracking, measuring 
momenta with a resolution of 5p±/p± of 0.4% at large angle and reconstruction of two 
track intersections, vertices, to an accuracy of ~3 mm. A superconducting coil around 
the EMC provides a 0.52 T magnetic field. The low-beta insertion quadrupoles are in the 
middle of KLOE. They are therefore surrounded by two compact tile calorimeters, QCAL 
Q, used as veto for otherwise undetected photons absorbed by the quadrupoles. 

The EMC is divided into a barrel and two endcaps covering 98% of the solid angle. 
Modules are read out at both ends by photomultipliers, PM, with a readout granularity 
of ~4.4x4.4 cni'^ for a total of 2440 cells. The calorimeter thickness is ~15 radiation 
lengths, Xq. Both amplitude and time information are obtained from the PMs. The signal 
amplitude measures the energy deposited in a cell and its time provides both the arrival 
time of particles and the position along the modules of the energy deposits, the latter by 
time difference. Cells close in time and space are grouped into a "calorimeter cluster". 
The cluster energy E is the sum of the cell energies. The cluster time T and position R 
are energy-weighted averages. R indicates the cluster position with respect to the detector 
origin of coordinates. Energy and time resolutions are ge/E = b.lVo/\/E (GeV) and 



at = 57 ps/yi*^ (GeV) © 100 ps, respectively. The photon detection efficiency is ~ 90% at 
E = 2{) MeV and reaches 100% above 70 MeV. 

The QCAL calorimeters, ~5Xo thick, have a polar angle coverage of 0.94 < jcos0| < 
0.99. Each calorimeter consists of 16 azimuthal sectors of lead and scintillator tiles. The 
readout is by wavelength shifter fibers and photomultipliers. The fiber arrangement allows 
the measurement of the longitudinal coordinate by time differences. 

Only calorimeter signals are used for the trigger ||l^. Two isolated energy deposits, 
> 50 MeV in the barrel and E > 150 MeV in the endcaps, are required. Identifica- 
tion and rejection of cosmic-ray events are also performed by the trigger hardware. A 
background rejection filter, Filfo |11], based on calorimeter information runs offline. Filfo 
rejects residual cosmic-ray, machine background and Bhabha events degraded by grazing 
the QCAL, before running event reconstruction. 



3. Search of Ks — >■ 77 with a pure Ks beam 
3.1 Ks tagging and event preselection 

The mean Ks and Kl decay lengths in KLOE are As ~ 0.6 cm and Al ~ 340 cm re- 
spectively. About 50% of the produced i^T^-mesons reach the calorimeter before decaying. 
-fC^-mesons are very cleanly tagged, with high efficiency ~30%, by identifying a Ki in- 
teraction in the EMC, which we call i^^^-crash. A i^^^-crash has a very distinctive EMC 
signature: a late {{(3k) — 0.22) high-energy cluster with no nearby track. The average value 
of the e^e~ collision center of mass energy, W, is obtained with an accuracy of 30 keV for 
each 100 nb~^ of integrated luminosity, by reconstructing large angle Bhabha scattering 
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events. The mean interaction point, IP, position and the p^f, momentum are also obtained. 
The value of W , and the JCi-crash cluster position provide, for each event, the trajec- 
tory of the Ks with an angular resolution of 1° and a momentum resolution better than 1 
MeV. In the analyzed sample, corresponding to an integrated luminosity J >Cdt=1.9 fb^^, 
we observe ~700 x 10^ tagged JC^-mesons. Using the most recent value of BR(iCs — > 77) 
Q, we expect ~1900 tagged Ks — > 77 events. Because of tagging, we have no Kl^2^ 
background, the major contamination in the NA48 measurement. The main background 
in our analysis is due to Ks 27r^ events with two photons undetected because out of 
geometrical acceptance or not reconstructed in the EMC. 

We estimate all backgrounds with the KLOE Monte Carlo, MC, ||T^. We produced 
(j) decays to all channels corresponding to an integrated luminosity J>Cdt~1.5 fb^^. In 
addition, for the signal we use a very large sample of MC Ks — > 77 events, equivalent to 
f Cdt'^100 fb~^. In the simulation, the photon detection efficiency and resolutions have 
been tuned with data using a large sample of tagged photons from <p —>' '7r"'"7r~7r'^ events 
selected using only drift chamber information |11|. Kl interactions in the EMC are also 
simulated. 

Since the Ks decay length is approximately 1/10 the distance traveled by a photon 
in our time resolution we take all -ftT^-decay photons as originating at the IP. A prompt 
photon is defined as a neutral cluster in the EMC, satisfying the condition |T — R/c\ < 
min(5<Tj, 2ns), where T is the time of flight (TOE) and R = |R| indicates the cluster 
position with respect to the detector origin of coordinates, at is the total time resolution. 
After tagging, we define a signal-enriched sample by requiring two and no more than two 
prompt photons in the event. While the minimum energy of photons from Ks^JJ is 
197 MeV, photons from i^5^27r°^47 can be much softer, > 15.8 MeV. Also at this 
momentum our resolution is of 0(40%). To maximize Ks^2tt^ rejection we therefore 
consider all clusters with E >7 MeV, and |cos(0)| < 0.93. The distribution of photons 
from Ks^2'7T^ not detected by the EMC is peaked at |cos0|=l, as shown by the MC 
spectrum in Eig. |^. 
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Figure 1: Angular distribution of photons from Ks^2tt'^ with two photons in the EMC. Recon- 
structed photons solid-line histogram, undetected photons points. 



After these cuts, we are left with 550,000 events, a signal efficiency of ~83% and a 
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signal over background ratio S/B~l/300. The background is mostly from Ks^'2^t^^ events 
(99.1%) and a 0.7% contamination of false Ki-cias\i from K^K~ events. There is also a 
residual background from i^^^-decays other than 0.2% from vr+vr", 0.02% from -nlu. To 
improve background rejection, we veto events with photons absorbed by the QCAL. Fig. ^ 
shows the distribution of the difference between the reconstructed and the expected time 
of the QCAL signals, ATq. The in-time peak is due to Ks^^tt^ with photons reaching 
the QCAL. 

The oscillating distribution is due to 
machine background events and shows 
the period of the beam bunches. All 
events having at least one hit in QCAL 
with energy above threshold and in a 
time window, TW, defined by |ATq| < 
5 ns are vetoed. This veto removes ~ 
70% of the background, while retaining 
high efficiency for the signal. The signal 
loss is ~0.04%. 

We must however correct for the sig- 
nal loss due to the accidental coincidence 
with machine background signals in the 
TW. The correction is Cq = 1 — Pq^ , where -Pq^ is the probability of a random coin- 
cidence in the TW. The latter is taken as the average of values obtained in two different 
out-of-time windows, one early and one late with respect to the collision time. We estimate 
the systematic error from the value of obtained from reconstructed iCs— >7r+7r~ decays 
where no photons are present. We find: P^^ = (3.51 ± 0.04stat i0.26syst)%- At the end of 
the acceptance and QCAL veto selection, we remain with 157 xlO^ events. The S/B ratio 
is ~l/80 at this stage. 

3.2 Kinematic fitting and event counting 

To improve the S/B ratio, we perform a kinematic fit imposing seven constraints: energy 
and momentum conservation, the kaon mass and the two photon velocities. Input variables 
to the fit are the IP coordinates, the Ks decay point, the Kg momentum |p|, the interaction 
points of the two photons in the EMC and the two cluster energies. All of these 15 variables 
are adjusted by the fit. There is no unmeasured variable to be determined. So this is a 
7-C fit with the number of degrees of freedom being dof=7. Fig. 3a, 3b and 6a, show a 
peak in at ~5 as expected for dof=7. 

Fig. ^ shows the \^ distribution from the fit for data and MC events, after acceptance 
selection, before and after applying the QCAL veto. The background has high values. 
Rejecting events with > 20 we retain ~ 63% of the signal while considerably reducing 
the background. The S/B ratio improves from 1/80 to 1/3. After this cut, the background 
is entirely due to Kg — > 2tt^ events with two undetected photons. Background, Fig. can 
be further reduced using the 77 invariant mass My^, and the photon opening angle in the 
kaon rest frame, 6* . Since the kinematic fit imposes the kaon mass as a constraint, we 
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Figure 2: Inclusive distribution of the difference 
between the measured arrival time and the expected 
time of flight of hits in QCAL for events tagged by 
a i^L-crash with two prompt photons (solid line) or 
with a reconstructed Ks — > tt+tt^ decay (points). 
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Figure 3: distributions for tagged Ks events witli two prompt pliotons: before (a) and after 
(b) QCAL veto. 
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use the measured variables values before fitting. Fig. ^ shows plots of M-y^ vs cos 6'^^ for 
data, MC background and MC signal events. 

To check the MC description of 
the EMC function of the pho- 

ton energy, we inspect the energy pulls 
of the kinematic fit for 2it^ Ks de- 
cays. We use a data sample corre- 
sponding to / £di~80 pb~^ and equal 
MC statistics. We select tagged Ks- 
mesons and ask for four prompt pho- 
tons. An energy scale correction of ^ 
1.02 is required to improve the match 
between MC simulation and data. Af- 
ter applying this correction, the MC 
ability to reproduce signal spectra is 
tested with a control sample of i^L— > 
77 events decaying near the beam pi- 
pe, with the i^T^-meson tagged by a 
well reconstructed ^C^— >7r+7r~ decay. 
The BR for Kl^2j is 5.74 x 10"^ 
which together with the lifetime, Tjy- 
=5.08 X 10^*^ corresponds to an equiv- 
alent BR(i^j^^27) of 1.6 X 10~^ per cm of path. Thus decays within 30 cm of the IP 
provide a sample of Kl^Ij larger than that of Ks^2'y and with a background level from 
27T^ decays smaller by three orders of magnitude. The vertex position is calculated by 
knowing the Kl flight direction and the time of flight of the two photons with a precision 
of ~ 1.5cm. 
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Figure 4: Scatter plot of M^^ vs cos 6**^, after pre- 
selection, for data (a), MC background (b) and MC 
signal (c). The solid curve represents the signal domi- 
nated region. 
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Data (MC) corresponding to fCdt=200 (450) pb~^ are used. Events are selected as 
for the Ks 27 decays, including the kinematic fit. The background is negligible after 
requiring < 20. A gaussian fit to the M-y^ distributions is shown in Fig. H. Data and 
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Figure 5: 77 invariant mass for the Kl 77decays near the beam pipe. 



MC energy scales agree to better than 0.2%, ~l/5 of the error which is quite satisfactory. 
The resolution agrees to 2%. The and cos 6**^ distributions of the events. Fig. ^ a 
and b, confirm the simulation results for Ks^2j decays. 




Figure 6: Distributions of (a) and 0086**^ (b) for Kl^jj decays near the IP. Black points are 
data, grey histogram is the MC simulation. The plot of has been done with preselection cuts on 
the two photons and a cut on cos 6**^ below -0.998. The cos6'*^ distribution required a cut at 
20. 

To obtain the number of A'5^27 events, we perform a 2 dimensional binned-maximum- 
likelihood of the the final sample distribution in the M^^ and cos 9*^ variables. The like- 
lihood function uses the MC generated signal and background shapes taking into account 
data and MC statistics. The fit gives N(77) = 711 ±35, with a xV^of = 854/826. The fit 
CL is 24.3%. 

Projections of data and fit are shown in Fig. ^. The signal cos 6*^ distribution is 
peaked at cos0=— 1 while the M^^ distribution is gaussian at the Ks mass. The back- 
ground is less peaked at cos0=— 1 and lower and broader in mass. As an independent 
check of the fit quality, we show in Fig. ^.a the distribution for data and MC after min- 
imization. A similar comparison is done also for the angular photon spectrum (Fig. §.b). 
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Figure 7: Distributions ofcos(6'* ) (a) and Af^^ (b) for the final sample. 




Figure 8: Distributions of (a) and inclusive cos0^ of the two photons in the event (b) for the 
final sample. 



which clearly indicates the presence of a flat component due to signal, as expected for the 
two body decay of a spin particle. 



4. Branching ratio evaluation and systematics 

The branching ratio is obtained from N{Ks — > 77) using for normalization the yield for 
Ks^2Tr^ in the same sample of tagged Xs-mesons by counting events with four prompt 
photons: 

■^■'(^^ - = nIk" X ^"°"f2'?,' X - 2.») X R. (4.1) 

N{Ks 27rU) eTOT(77|tag) 

The total efficiencies have been evaluated by MC after Ki-ciash. tag. The signal total 
efficiency is the product of the efficiencies for the acceptance selection, the QCAL cut and 
the cut: 

eTOT(77) = esei(77) x eQ(77) x (^x^ill)- (4-2) 
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For the normalization sample, the efficiency is related only to the acceptance of four 
photons. The ratio, R^, of all other efficiencies (triggering, Filfo filter and tagging) be- 
tween signal and normalization sample should be identically one. From MC we find 
i?e = 1.001 lb O.OOlstat- The difference from unity is added as contribution to the sys- 
tematic error on the BR. 

For the signal selection efficiency we find: 

eseiill) = (82.9 ± 0.2stat ± 0.2syst)%. (4.3) 

The large selection efficiency is due to the wide angular coverage of the calorimeter, the 
low energy threshold used and the almost flat angular distribution of the decay products. 
The systematic error assigned to this efficiency has been found by varying the data-MC 
correction of the cluster reconstruction efficiency. The efficiency for the QCAL cut is 
found from MC to be Eq''^ {2'y) ~ 99.96%. Applying the correction due to accidental losses 
described in sec. ^T] we obtain: 

eQ(27) = e^C(27) x Cq = (96.45 ± 0.04,tat ± 0.26syst)%. (4.4) 

The MC efficiency of the cut is e^2 = (63.3 ± 0.7)%. The systematic error related to 
the knowledge of the data— MC difference in the scale has been evaluated by using the 
Kl 77 control sample. For the chosen cut, we evaluate the data over MC ratio, R^, 
of the x^ cumulative distributions and we get {R^ — 1) = (— 0.5± 1.8)%. We conservatively 
assign the error on R as the contribution of the x^ scale to the systematic error. 



Source 


+ABR/BR (%) 


-ABR/BR (%) 


Trigger, Filter, Tag 


0.10 


0.10 


Signal acceptance 


0.17 


0.17 


QCAL veto 


0.02 


0.26 


X^ scale 


1.80 


1.80 


Background shape 


1.04 


0.98 


QCAL TW change 


0.53 


0.49 


X^ change 


0.99 




MC Energy scale 




0.79 


2D-Fit binning 


0.96 


0.98 


Normalization sample 


0.15 


0.15 


Total 


2.56 


2.48 



Table 1: Breakdown of the contributions to the total systematic error for the BR(if5 — > 77). 

The systematic uncertainties connected to the signal counting have been evaluated by 
repeating the analysis and the fit in different ways. The most delicate point is related to 
the simulation of the background shape. The MC shows a good agreement with data for 
background-enriched samples obtained by requiring a complementary cut on x^j such as 
30 < x^ < 500. Moreover, to test the fit stability in different regions of the M^^ - cos 9*^ 
plane, we have determined how much the result varies when: (1) reducing the fit-region 
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along the cos^*^ axis moving the lower boundaries from 0.999 to 0.9995 or (2) fitting only 
in a signal dominated region shown by the ellipse in Fig. ^. The maximum variation of the 
BR for these tests is reported as background shape in Tab. ||. 

We have also tested the stability of the branching ratio when modifying the width of 
the time window used for the QCAL veto from ±5 ns to ±4, ±6 ns. Similarly, the cut in 
'X^ has been changed from 20 to 10 and 24. We have then repeated the fit by applying to 
the MC an energy-scale correction of +0.4%, a factor of two larger than what measured 
with the Ki^ — > 77 control sample. We have also checked that regrouping the bins of the 
2-D plot by factors from 2 to 5 does not modify substantially the result. For all of these 
cases, the maximum variation of the BR obtained is used as systematic error and shown 
m Tab. |. The sum in quadrature of all entries is used as total systematic error. 

For the normalization we count Ks 27r^ tagged events with four prompt photons. 
An efficiency of 

esei(27r0) = (65.0 ± 0.2stat ± 0.1syst)% (4.5) 

is found by MC. As for the signal, the systematic uncertainty related to the cluster de- 
tection efficiency is evaluated by varying the data-MC correction curves. After correcting 
for esei(27r°), a number of (190.5 ± 0.2) x 10^ Ks 27r'^ tagged events is obtained. The 
systematic uncertainty related to the presence of machine background clusters, fragmen- 
tation and merging of clusters is estimated by repeating the measurement in an inclusive 
way and counting tagged events with 3, 4 and 5 photons. The overall systematic error for 
the normalization sample is reported in Tab. [l|. 

To evaluate BR{Ks 77) we use the latest PDG |l2| value BR{Ks 27r^)= (30.69 ± 
0.05)%. See also |l|]. We obtain: 

BR{Ks ^ 77) = (2.26 ± 0.12,tat ± 0.06syst) x 10^^ (4.6) 

We have repeated the measurement by subdividing the data in two sets to check stability 
for the slightly different running conditions: 1) 0.4 fb"^ from 2001-2002 and 2) 1.5 fh~^ 
for 2004-2005. Also the simulation has been divided accordingly. We get BR{Ks 77) = 
(2.24 ± 0.30stat) X 10"*^ in 2001-2002 and BR{Ks 77) = (2.26 ± 0.13stat) x 10"^ in 
2004-2005, which are in excellent agreement. 

Fig. m shows our result and other existing measurements of BR(i^5— >77) as well as the 
0{p'^) theoretical prediction. There is a 3 u's discrepancy between the present result 
and the measurement of NA48. 



5. Conclusion 

From ~ 2 billion </> mesons collected with KLOE at DA<I>NE, we have measured the 
BR{Ks 77) with a 5.3% statistical uncertainty and a ~ 2% systematic error. We 
obtain a BR result which deviates by 3 (t's from the previous best determination. Precise 
xPT theory calculation for this decay are done at 0{p^). Higher order effects are predicted 
to be at most of the order of ~ 20% of the 0{p'^) decay amplitude. Our measurement is 
consistent with negligible higher order corrections. 
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measurements and xPT predictions. 
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